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Abstract 

In this paper, we propose a generic and systematic approach for 
study of the electronic structure for atoms or molecules. In partic- 
ular, we address the issue of single particle states, or orbitals, which 
should be one of the most important aspects of a quantum many-body 
theory. We argue that the single-particle Green function provides a 
most general scheme for generating these single particle states or or- 
bitals. We call them the correlated atomic or molecular orbitals to 
make a distinction from those determined from Hartree — Fock equa- 
tion. We present the calculation of the single particle properties (i.e., 
the electron affinities (EA's) and ionization potentials (IP's)) for the 
H2O molecule using the correlated molecular orbitals in the context of 
quantum chemistry with a second-order self energy. We also calculate 
the total ground state energy with a single Slater wavefunction deter- 
mined only from the hole states. Comparisons are made with available 
experimental data as well as with those from the Hartree — Fock or 
density functional theory (DFT) calculations. We conclude that the 
correlated atomic or molecular orbital approach provides a strictest 
and most powerful method for studying the single-particle properties 
of atoms or molecules. It also gives a better total energy than do the 
Hartree — Fock and DFT even at the single Slater determinant level. 



1 



It promises that a correlation theory based on the correlated atomic 
or molecular orbitals will become an approach which possesses the ad- 
vantages and also overcomes their shortcomings of current quantum 
chemistry methods based on either the conventional quantum many- 
body theory or the DFT. 

1 Introduction 

The single particle approximation, or the concept of atomic or molecular 
orbitals in the context of quantum chemistry, is a natural and almost a nec- 
essary scenario for solving an interacting many-electron system for atoms, 
molecules, or solids [T]. This is a reflection of not only a physical existence 
but also possibly a mathematical reality. The usual equation which is being 
used to determine the orbitals is the Hartree — Fock equation j21 IS]- The 
rest of endeavor to remedy the approximation resulting from a replacement 
of the whole many-body wavefunction by a single determinant used in the 
HF scheme is called the correlation issue. This is a most difficult problem 
and constitutes the major activity of researches for the quantum chemistry 
community in the last 50 years 4J . According to the energy scale principle we 
described in paper [5], the Hartree — Fock scheme should be a good approxi- 
mation when the single determinant wavefunction dominates and there is no 
any significant mixing with the nearby configurations. This is typically the 
case when we compute the energetics for molecules with a stable geometric 
structure. The subsequent perturbation correction for the correlation such 
as MPPT is also proved to be powerful However, there are the situations 
when the configuration mixing is a prominent or dominant phenomenon, and 
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the description with more than one single configuration seems necessary. This 
includes the calculation of transition states or excited states, and for open- 
shell molecules, etc. The computation based on the Hartree — Fock equation 
has shown to be insufficient, and the corresponding perturbation correction 
has proved not to be convergent [7| |H1 E]- The MCSCF approaches have 
been introduced to investigate this type of nondynamic or static correlation 
issue, and they have become one of the most popular approaches for the 
study of molecular electron correlation 0021 El E|- However, the size of 
molecular systems that this type of approaches can address are still limited 
because of the difficulties in selecting the appropriate configuration states 
and in achieving the convergence to the correct state of the interest fT^ . 
Another important and significant advance in the fields of electron correla- 
tion is the development of density function theory (DFT) [t^ I14j. Instead 
of working with a multi-configurational framework, it intents to incorporate 
the exchange-correlation effect into a single-particle potential formalism. It 
has already shown its very usefulness in the study of the electronic structure 
for large systems with utilization of relatively smaller computational efforts. 
Nevertheless, there exist some serious drawbacks for the method when seen 
either from theoretical consideration or from the practical performance in 
calculation. One shortcoming is that the theory can only study the ground 
state problem, and cannot treat the same eigenstate problem for excited 
states within one theoretical framework. Another serious problem is that the 
actual form for the exchange- correlation is unknown, or the theory itself gives 
no clue for how to approach it. Moreover, the approach fails to or can not do 
the accurate computation for the points or situations when the configuration 
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mixing is important [1^1 CHI HI] • Indeed, it should be a very difficult thing, 
intending to replace the intrinsic many-body effects such as static correlation 
or configuration mixing with a sm^^/e-particle formalism. 

Recently, we have demonstrated that a general quantum many-body per- 
turbation theory can not only be used for understanding the various elec- 
tronic phenomena including the nature of chemical bonds but also serve as a 
unified theme for constructing general electronic structure theories and cal- 
culation schemes. This also includes the study of important issues of electron 
correlation [B]. This pinpoints the direction and paves the way for the fu- 
ture investigation. In this paper, we add another important ingredient to 
the field of electron correlation or electronic structure theory in general. We 
emphasize our investigation on the issue of single particle approximation, 
or the atomic or molecular orbitals for the quantum chemistry calculation. 
From the perturbation point of view, this corresponds to defining a refer- 
ence Hamiltonian ^] . We will show that there exists a strict theoretical 
formalism, called the single-particle Green function, which provides a most 
general scheme for generating or determining these single-particle states up 
to present time. The theory of single-particle Green function has been de- 
veloped for a long time and used in many different ways but its full physical 
meaning or context is not totally understood or appreciated. This paper aims 
at a beginning for a systematic investigation of electron correlation based on 
the single-particle Green function formalism and within the quantum many- 
body perturbation theory [5^. In the next Section, we present its definitions 
and equations in both time and energy domains. In particular, we give an 
energy eigenequation that solves the single-particle states. We analyze its 
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intrinsic structure and compare it with other methods. In Section 3, we cal- 
culate both the single-particle properties and the total energies for the H2O 
molecule using the Hartree — Fock, DFT, and correlated molecular orbital 
approaches. In the final Section, we analyze and discuss our results for the 
calculations and also do the comparison with each other including the corre- 
sponding experimental data. We also propose a generic electronic structure 
theory and outline the future research. 



2 Theory 

Two time {t,t') and single-particle (or hole) Green function is defined as fHl 
120] 

Gixt,xt) = -^(^o|T{^(x,t)^+(x ,t')}|^o), (1) 

where T is Wick time-ordering operator, and ip{x,t) and '?/'^(af', t') are the 
field operators in the Heisenberg picture associated with the coordinates x, 
which includes both spatial r and spin x degrees of freedom. The |\E'o) is the 
exact ground state of an A^-electron system being studied. Its Hamiltonian 
in the field operator representation can be written as 

H = j 'ip^{x)h{x)'tp{x)dx + j 'ip^{x)'ip^ (2) 

where the one-body operator h{x) is the sum of the electronic kinetic energy 
operator and its interaction with the nucleus 

m = -T^V^ - J2 Zpv{x, Rp), (3) 
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and the two-body operator v{f, r ) is the Coulomb potential 



v{f, r ) 



1 



(4) 



r — r 



In the energy domain, the Green function takes the form 




(5) 



where 



Mx) = {^omx)\^n{N+l)), er, = Er,{N + l)-Eo for en > II, (6) 



= (*„(7V - l)|V'(f)|*o), en^Eo- En{N - 1) far e„ < (7) 



The wavefunctions |^„(A^±1)) and energy levels En[N ±1) are for the A^±l 
electronic systems. The functions {0„(x)} are the ones of sm^'Ze-particle 
coordinates, and are called the particle states for those defined by Eq. (6) 
(cn > yw), and the hole states for those defined by Eq. (7) (e„ < /x), where 
/X is the chemical potential. The corresponding energy e„ are the electron 
affinity or the electron ionization potential, respectively. A very important 
feature of these single-particle states is that they form a complete 

set as shown below. 



where n is for all the hole or particle states. The Eq. (5) is called the 
Lehmann representation. 

Define the average classical Coulomb potential by 



or 



(8) 



n 




(9) 
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where 

p(f) = (^o|V^+(x)^(x)|*o), (10) 

is the one-electron probabihty density, then the Green function in the energy 
domain satisfies the following equation, 

{e — h{x) — V{x)} G{x, X ] e) — J Tj{x, X ] e)G{x , X ] e)dx =6{x — x), (11) 

where the operator T,(x,x ; e) is called the self-energy operator which is non- 
local and energy dependent. From this equation for the single-particle Green 
function and its Lehmann representation (5), we can get an equation that 
the single-particle states {0„(a;)} satisfy 

{h{x) + V{x)} 0„(f) + J S(x, x; en)(t)n{x)dx = e„0„(x), (12) 

or 

{/i + V + S(e„)}|0„) =e„|0„) (13) 

in a more general Dirac notation. It is called the Dyson equation or the 
energy eigenequation for the quasi-particles in the current literature [201 
1221 12n] • When we do the comparison with the Hartree — Fock equation or 
the Kohn — Sham equation [21 El US] , it seems that the self-energy operator 
S is related to the exchange and correlation effects of an interacting many- 
electron system beyond that of the classical Coulomb interaction. Unlike 
the Kohn — Sham equation, however, where the explicit analytical potential 
for the exchange-correlation potential is unknown, the self-energy operator 
has intrinsic structure, and, for example, can be expanded as a perturbation 
series as follows, 

S = S(°) + sW + ... + S(") + .... (14) 
7 



They have exphcit physical interpretations and therefore can be approached 
in a systematic way[2niinil22ll23l21l2Sll2nil2Zll21l2i- Another important 
feature of Eq. (12) is that the single-particle states are defined for both hole 
state (Eq. (6)) and particle states (Eq. (7)), and therefore there exists the 
concept of a fundamental excitation in the present formalism. In other words, 
we can form the configurations based on these single particle states. Further- 
more, since they constitute a complete set of single-particle states, as shown 
in Eq.(8), any A^-electron wavefunctions can be expanded as a linear combi- 
nation of these configurations. For these reasons, we can regard the equation 
(12) as a most general eigenequation for creating the single-particle states 
or the atomic or molecular orbitals at present time. It is the corresponding 
one-particle description of an interacting many-hodj system For 
clearness and easiness to be understood, we call the single-particle states 
determined by Eq. (12) as the correlated atomic or molecular orbitals in 
order to make a distinction from those determined from the Hartree — Fock 
equation. Obviously, they will catch the full Hamiltonian (2) more than do 
the Hartree — Fock orbitals. 

The successfulness for obtaining the most appropriate correlated atomic 
or molecular orbitals {0„(a;)} will depend on how well we can obtain the 
correct self-energy operator S. This will in turn depend on what kind of 
wavefunctions or what level of theories we select as the reference or the ini- 
tial wavefunction for our construction of S since the Eq. (12) is an iterative 
equation for determination of {e„} and {0„(s)}. Obviously, there will be dif- 
ferent choices for different species or for different molecular geometries being 
studied as have already been demonstrated in many existing quantum chem- 
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istry calculations. Several types of perturbation schemes for the self-energy 
operator have already been developed either from solid state physics com- 
munity or by quantum chemists |2niinil22ll23l23l23l2Hll2Zll2Hll23- These 
include the functional derivative method [201 122 12^] , the superoperator 
formalism f2M |2S1 12S1 12Z] , the diagrammatic expansion method j2H] , and the 
equation of motion approach [21]. 

3 Calculation and Results 

In this section, we present the computation of the single-particle properties 
and total energies for H2O molecule. We employ the Hartree — Fock method, 
DFT, and correlated molecular orbital approach we describe above for the 
calculation and do the corresponding comparison. 

The geometric parameters for the water molecule are taken from experi- 
mental observation which are i?(0-i^) = 0.957i, and IHOH = 10A.5{deg) 
For the Hartree — Fock calculation, we use the cc — p VTZ basis set j22] ■ The 
calculated energies for the first ten molecular orbitals are listed in the second 
column of Table 1. The computed total energy is shown in the Table 2. For 
the DFT calculation, we use the same set of basis functions. The exchange- 
correlation functional is approximated with the B3LYP scheme [33 Ell- The 
result for the first ten Kohn — Sham orbital energies is listed in the third 
column of the Table 1. The total energy is shown in the Table 2. For the 
computation based on the correlated molecular orbitals, we take the second- 



9 



order approximation for the self-energy operator, 

S,,(E) = E«(E) + sg)(E). (15) 

The detailed forms for the self-energy operator with different orders are 
dependent upon the reference states chosen ESI- For the closed-shell 
molecules, if we pick the Hartree — Fock orbitals as the reference states, 
the first-order self energy vanishes, and the second-order self-energy is given 

by UHl 

Y-(2)/77^ _ {rs\ia) {2{ja\rs) - {aj\rs)) ^ {ab\ir) {2{jr\ab) - {rj\ab)) 

^ij y^) - 2^ ^r— — - — + TT— — - — 

ars -C/ + Ca — Cr — abr ^ + ^r — — 

(16) 

where a,b,... are the spatial hole states, and r, s, ... are the spatial parti- 
cle states. If we choose the Kohn — Sham orbitals as the reference states, 
however, the first-order self-energy takes the form 

N 

sgH^) = -(^|K.c|j)-EHaj), (17) 

a 

and the second-order self-energy remains the same as that for the case of 
the Hartree — Fock orbitals. We solve the eigenequation (12) for the quasi- 
particles with the cc — pVTZ basis set. When the Hartree — Fock orbitals 
are used as the reference state, the calculated quasienergies for the first ten 
correlated molecular orbitals are shown in the third column of the Table 1. 
The resulting total energy with the single determinant using the first five 
doubly-occupied hole states is also listed in the table 2. When the DFT 
determinant is employed as the reference state, the corresponding results are 
listed in the forth column of table 1 or table 2. All the computations are 
done with the Hondo — v99.6 suite 
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4 Discussion and Conclusions 



In this paper, we present a novel approach for the study of electronic structure 
of atoms and molecules related to the single-particle Green function theory. 
We argue that the single-particle Green function provides a most general 
theoretical framework for generating the atomic or molecular orbitals for the 
atoms and molecules. Based on this statement, we have calculated both 
the energies of these single-particle states and total energies for the H2O 
molecule [SZlEHlEn]- For the total energy, a single-determinant wavefunction 
composed of hole states only is used for the computation. At the same 
time, the calculations are also performed with the Hartree — Fock and DFT 
methods. 

When compared with the experimental ionization energy or electron affin- 
ity for H2O molecule we see that the correlated molecular orbitals with 
the Hartree — Fock orbitals as the reference state gives the better results 
than the ones from the Hartree — Fock or DFT methods. The total energies 
obtained with three different methods are also compared to the one obtained 
from the experimental observation pT|l37|l3Sl l39j. The correlated molecular 
orbital approach results in the best value. Of course, the calculation can be 
further improved by choosing the DFT as a reference wavefunction. We have 
the similar conclusion. 

Since the work of Heitler and London in the calculation of the electronic 
structure for H2 molecule, which is the indication of the beginning of the field 
of quantum chemistry, it has the history of development for more than eighty 
years. However, there is a fundamental issue, i.e., the quality of atomic or 
molecular orbitals, which has been neglected for a long time. This paper 
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addresses this "quality" issue for single-particle states or orbitals in many- 
body theory. From the perturbation theory point of view, this corresponds to 
a definition of the reference Hamiltonian, which is crucial in the minimization 
of dynamic correlation energy or convergence of perturbation series. It is 
also critical in providing the best single particle properties. Both of the 
calculated single-particle properties and total energies have explicit physical 
interpretation and are subject to the test from experimental observations ^T] . 

From above analysis, it is obvious that when the concept of correlated 
atom or molecular orbital is incorporated into the quantum many-body per- 
turbation or coupled cluster theory, it will provide a most powerful quan- 
tum many-body approach for the study of electronic structure of atoms or 
molecules. On one hand, its single-particle properties have obvious physi- 
cal meanings which is in contrast to the case for the DFT. Furthermore, it 
can go beyond the single-determinant level and form configurations. There- 
fore, it can study the issues when configuration mixing is important. On the 
other hand, when doing the comparison to the traditional quantum many- 
body theory based on the Hartree — Fock or MCSCF orbitals, the correlated 
orbital method not only has provided a better description of single-particle 
properties, but also gives us the better convergence at the configuration level 
and therefore provides a more powerful computational scheme. For these 
reasons, we could claim that the correlated atomic or molecular approach 
will be a most general ah initio correlation method for electron structure cal- 
culations. It possesses the advantages and also overcomes their shortcomings 
of current DFT and conventional correlation approaches based on the atomic 
or molecular orbitals determined from the Hartree — Fock or MCSCF . 
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Of course, it has been a very difficult task for a long time to get the ap- 
proximate self-energy operator to the higher orders. However, the intrinsic 
structure such as its perturbation series expansion has offered us a possi- 
bility instead of an outside model for the approximation. Furthermore, the 
further study of this underlying intrinsic structure will tell us more universal 
things which might be true even for a many-body theory or system in gen- 
eral. Henceforth, the continuing investigation of the higher order self-energy 
operators and their relations will be a rewarding research |^ . 

An interesting point needed to be mentioned is that the self-energy oper- 
ator in Eq. (12) does not have to be Hermitian which corresponds to the sit- 
uation when ipn{x) is a real orbital. Here the imaginary case for the operator 
is related to the electron dynamics which is left as a future investigation [201 

Finally, if we fully explore the usefulness of the pseudopotential theory, 
combined QM / MM approach, or linear scaling algorithms and so forth, the 
correlation theory based on the correlated atomic or molecular orbitals will 
provide to us a most robust approach for the study of electronic structure 
even for large systems 
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Table Caption 



Table 1. The single-particle properties or orbital energies (in a.u.) of 
H2O molecule from the calculations based on the Hartree — Fock, DFT and 
correlated molecular orbital approaches as well as from the experimental 
measurement. 

Table 2. The total energies (in a.u.) of H2O molecule from the calcula- 
tions based on the Hartree — Fock, DFT and correlated molecular orbital 
approaches as well as from the experimental measurement or CI calculation. 
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